STABILITY IN RANDOM BOOLEAN CELLULAR 
AUTOMATA ON THE INTEGER LATTICE 



F. MICHEL DEKKING, LEONARD VAN DRIEL AND ANNE FEY 



Abstract. We consider one-dimensional random boolean cellular au- 
tomata, i.e., the cells are identified with the integers from I to N. The 
behavior of the automaton is mainly determined by the support of the 
random variable that selects one of the sixteen possible Boolean rules, 
independently for each cell. A cell is said to stabilize if it will not change 
its state anymore after some time. We classify the one- dimensional ran- 
dom boolean automata according to the positivity of their probability 
of stabilization. Here is an example of a consequence of our results: if 
the support contains at least 5 rules, then asymptotically as N — > oo 
the probability of stabilization is positive, whereas there exist random 
boolean cellular automata with 4 rules in their support for which this 
probability tends to 0. 



1. Introduction 

We study the dynamics of cellular automata as amply discussed in Stephen 
Wolfram's book "A New Kind of Science" (|12j). The cells, which are iden- 
tified with the integers from 1 to N, can be in two possible states or 1, 
and the state at the next time instant is a function of the states of the two 
neighbouring cells. For this function, also called (evolution) rule, there are 
16 different possibilities, see also Table HJ Although there is a chapter in [12] 
on random initial conditions for CA's, there is hardly anything on random 
rules, or even on the possibility to assign different rules to different cells (so 
called inhomogeneous CA's). In this paper we will give a classification of 
one-dimensional Random Boolean CA's. 

Our results are also of interest in the context of Random Boolean Networks 
introduced by Stuart Kauffman in 1969 ([3]). Here not only the rules and the 
initial condition are chosen randomly, but also the neighbours of a cell, i.e., 
the N cells obtain a random graph structure. For each cell i independently 
a mother Mil) and a father F(i) are chosen uniformly from the remaining 
pairs of cells. This gives a graph with vertices i and directed edges (M(i),i), 
(F(i),i) for i = 1,...,N. Next, each cell % is assigned independently a 
Boolean rule $j uniformly from the set of 16 Boolean rules (i.e., maps from 
{0, l} 2 to {0, 1}, see also Tabled]). To start, all cells obtain independently 
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an initial value Xi(0) from the set {0, 1} with equal probability. The state 
of a RBN at time t will be denoted by X(t). Then the dynamics is given by 
the map 

X(t + l) = **(X(t)) 

defined by 

Xi(t + 1) = $i(X m {t),X Fii) (t)) 

for all cells i. 

A 1-dimensional random Boolean cellular automaton can be considered as 
an RBN with 

M(i) =i-l and F(i) =i + l, 

fori = 2,..., JV-1. Weput M(1)=JV,F(1) = 2, and M(A) = A-l, F(A) = 1. 
In Figured] we display a realization of the evolution of a RBCA in the usual 
way: zero's and ones are coded by white and black squares, space is in the 
horizontal direction, time in the vertical direction. 




FIGURE 1 . Evolution of a one-dimensional RBCA. After two time 
steps, most cells already become periodic. After five steps, the whole 
CA is periodic. 



We say a cell stabilizes if its value is not changed anymore after some time. 
In Figure [1] the stable cells yield the black or white stripes; in this RBCA 
example about 56% of the cells stabilize. 

Kauffman ([!]) found in his simulations that about 60-80% of the cells in 
a two-connected RBN stabilize. It was later shown rigorously by Luczak 
and Cohen ([5]) that in the limit as the number of cells N tends to infinity 
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100% of the cells will stabilize — size N = 10 000 in Kauffmans simulations 
is much too small to unearth this behavior. 

For the RBN's there is no natural model for the infinite size (N — > oo) 
network. For a one- dimensional CA there is: the graph is simply the set of 
all integers with edges between nearest neighbors. It therefore makes sense 
to speak of the probability cr* of stabilization of the infinite automaton. The 
result of Luczak and Cohen leads to the question whether this a* = 1. In 
this paper we will show that this is rarely the case, and moreover that there 
are many cases where actually a* = 0. By different cases we mean different 
distributions on the set of 16 Boolean rules. We remark that the results in 
[I] and [5] are restricted to the uniform distribution on the set of rules — see 
[21 El O [8] for some results for RBN's with other rule distributions. 
The second problem we will discuss is whether — > when N — > oo. We 
will show that this will often be the case, but may also fail to hold. 
For some previous work with rather partial results on 2-dimensional versions 
of RBCA's see e.g., pp. We will give an application of our results at the end 
of Section H] to a recent paper on RBN's of Greil and Drossel (0). 
Finally we mention that our interest in RBCA was also motivated by certain 
CA models of the subsurface in the geosciences (see [TTj). 



2. Random Boolean cellular automata 

Before we define RBCA, we describe their realizations which are inhomoge- 
neous cellular automata. The dynamics in such a CA is defined by attaching 
a Boolean rule to each cell. The rule is given by specifying what its value 
is on the four different values (1, 1), (1, 0), (0, 1) and (0, 0) of the neighbors 
of a cell. In Table (H all 16 different rules are defined. 

For ease of notation these are indexed as <pj where j is the number in binary 

notation obtained from the four bits in the row describing ifj . Anytime we 

refer to one of these 16 rules we use index j (note the different font); when 

we use index i we refer abusively to the rule assigned to cell %. 

An inhomogeneous CA's mapping ^ : {0, 1}^ -> {0, 1}^ is determined by 

a vector of rules, if = (y?i, . . . ,(p N ). Here for all time steps, <pi is the rule of 

cell i. Hence given a state x — (x\ . . .Xn) we get (^<p(x))i = (Pi{%i-i, 

for i — 2, . . . , N— 1. To deal with the problem that occurs in the end points 

we consider the automata as defined on a circle; a convenient way to express 

this is to define Xi + kN — x% for i = 1, . . . , N, and all integers k. 

To address ^ we will also write the vector of lower indices of the </?j (1 < 

i < N) mentioned. E.g., mapping (0,3,3,15) describes ^ : {0, l} 4 — > {0, l} 4 

with rule vector (ip , ip 3} <p 3 , ip 15 ). 
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Table 1. Definition of the 16 Boolean rules. The last col- 
umn gives a formula for </?j (x, y) when ip^ is afline. 



j 


^(0,0) ^(0,1) 


^(1,0) 


y?j(M) 


afline 




















1 


1 











no 


2 





1 








no 


3 


1 


1 








x + 1 


4 








1 





no 


5 


1 





1 





y + 1 


6 





1 


1 





x + y 


7 


1 


1 


1 





no 


8 











1 


no 


9 


1 








1 


x + y + 1 


10 





1 





1 


y 


11 


1 


1 





1 


no 


12 








1 


1 


X 


13 


1 





1 


1 


no 


14 





1 


1 


1 


no 


15 


1 


1 


1 


1 


1 



In a RBCA the rules tfi are chosen according to a distribution $ on the 
set of rules, independently for each %. The initial state is an i.i.d. sequence 
(Xi(0)) from a Bernoulli(l/2) distribution, independent from the Then 
the value of the i th cell at time t + 1 is given by 

(1) X i (t + l) = $ i (X i - 1 (t),X i+1 (t)). 

We denote by py the probability measure generated by the (Xj(0)) and 
the (<&j) on the space of infinite periodic sequences with period N. It is 
natural to consider also the probability measure P* on {0, 1} Z generated by 
two infinite independent i.i.d. sequences (Xj(0)) and (<3>j). 

3. Stability 

We say cell % stabilizes at time T if 

X i (t)=X i (T) foralU>T, 

and T is the smallest integer with this property. Let <tn be the probability 
that a cell in a random cellular automaton of size N with rule distribution 
<3> stabilizes. Because of stationarity, 



a N = (cell stabilizes) . 
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We first give an example were <jn does not converge as iV — > oo. 
Consider the RBCA's for N = 1,2... with rule distribution given by 

= <fe) = 1. 

Note that randomness is only involved because of the randomness in the 
initial state X(0) = (Xj(0)), and that ip e is the well known linear rule which 
can be written as (pe(x, y) — x + y. This addition, as all additions with O's 
and l's, is modulo 2. 

Proposition 1. For the ip 6 -RBCA one has 

liminf o"at = 0, lim sup = 1- 

Proof: It is well known (see e.g. [10J ) and easy to prove that the (p 6 - 
automaton satisfies 

Xi(t) = X^ t (0) + X i+t (0) if t = 2 P for some p. 

It follows that on = 1 for each N = 2 P which is a power of 2: 

X t (N) = X_ N+i (0) + X N+i (0) = 2X { {0) = for all i, 

which implies that Xi(t) = for all t > N. This proves that lim sup = 1- 
To obtain the statement on liminf, consider the subsequence of N of the 
form iV = 2 P + 1. Then we find for all % 

X t (N - 1) = X,(2 p ) = X_ 2P+i (0) + X 2P+i (0) = X i+ i(0) + X^O). 

From this we directly obtain that 

X (N) = X_ X (JV -1)+X 1 (N-1)= X_ 2 (0) + X (0) + X (0) + X 2 (0), 

and so 

(2) X (iV)=X_ 2 (0)+X 2 (0). 
Thus 

F N (X (N - 1) ^ X (N)) = F N (X-i(0) + *i(0) ^ X_ 2 (0) + X 2 (0)) = ^ 

Note that this computation does not only hold for t = N — 1, but for any t 
which satisfies t = —1 mod N. An infinite sequence of such t's is given by 

t = 2( 2fc+1 ) p , k = 0, 1, 2, For N of the form 2 P + 1 we therefore obtain 

that 

1 - a N = F N (X (t) ^ X (t + 1) for infinitely many t)>^. 

This inequality is a warm up for a more general result: we can generalize 
Equation (|2j) to for all m 

X (N + 2 m -2) = X_ 2m (0)+X 2m (0). 
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Since for m — 0, 1, . . . , £ the X„ 2 m (0) + X 2 m(0) are independent Bernoulli 
variables with success probability 1/2, it follows that 

F N (3m, 0<m<£: X (N + 2 m - 2) ^ X (N - 1)) = 1 - 1 

and employing the same infinite sequence of fs as for the £ = 1 case, we can 
deduce that liminfAr^oo = 0. □ 

It is convenient to define random variables Zj by Zj = 1 if cell « stabilizes, 
and otherwise. We then have (independently of i): 

<tjv = ¥ N (Zi = 1) , and a, = P*(Z, = 1) . 

The uniform RBCA is given by the rule distribution Pjv($ = = 1/16 for 
j — 0,1, ... , 15. 

Proposition 2. Let er^r fre the probability of stabilization of a size N uniform 
RBCA, a* that of an infinite size uniform RBCA. Then 

lim o~ N = a* and a* > 0. 

N^oo 

Proof: The two rules ip and ip 16 are called walls. We denote W = {ip ,ip 15 }. 
Cells that have obtained a rule from W are already stable at time or 1. 
We define for positive integers k and / the chambers Ck,i by 

C k , t = {<£ _ fc G W, £ W, $_ fc+2 £ W, . . . , £ W, $i G ^}. 

In addition we put C 0j o = {$o G Vr}. Since the chambers are disjoint events 
whose union has probability 1 we can write 

oo 

<7„ = P„(Z„ = 1) = P*(Z = 1 | C7 ,o) P*(C ,o) + ^ P*(Z = 1 I C fclI ) P *(Cm) • 

fe,«=i 

We will make a similar splitting for finite RBCA of size N. In case iV = 2M+ 
1 we consider the indices of the chambers from the set {— M, . . . , 0, . . . , M}, 
in case N = 2M from {— M+l, . . . , 0, . . . , M}. In the sequel we will assume 
N is odd, as our arguments can be transferred trivially to the case N even. 
Let D be the event 

M 

D= (C7 , U |J C h tf. 
k,i=i 

Then D occurs iff there is no wall among $_m, . . . , $o or among $ , • • • > ®m- 
Since Pv($; W) = P*($; £ W) = 7/8, we obtain that 

/7\M+1 /7n2M+1 /7\N/2 

W) = W) = Kl) -is) £2 (s) ' 

The crux of the proof is that for all 1 < k, I < M and for k = I = we have 
Pat(Z = 1 1 Ck,i) = P*(Z = 1 1 Ck,i) , 
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simply because both the finite automaton and the infinite automaton evolve 
between the walls independently of the evolution outside the walls. Thus 

\(Tn-v*\ < W*(Zo = l\D)F*(D)+F N (Z = l\D)F N (D) 



It follows that o n — > o~* exponentially fast. Note that 

on > P W ($ e W) = \, 

o 

and so <r* > 1/8 > 0. □ 

Actually we can obtain a good estimate of <t* using the exponential conver- 
gence of ctjv, by sampling uniform RBCA's for some large N and computing 
a 95% confidence interval. We found that ct 100 = 0.678058 ± 0.000029. In- 
terestingly, an exhaustive enumeration for iV = 6 yields already a value oe 
close to 0.68. 

We end this section with the remark that the 'wall'-property directly implies 
that the sequence (Zj) is strongly mixing under P*: 

P*(Z Q = £ Q , Z t = £l ) -> P*(Z = £ Q ) F*(Z = £l ) as i -»• oo. 

Therefore the ergodic theorem applies, yielding that we also have pathwise 
about 68% of stable cells. 

4. Stability by means of impermeable blocks 

Note that the proof of Proposition [2] will still be valid if $ is not uniformly 
distributed, but has an arbitrary distribution with Pat(<3> = ipj) > for all j 
(of course 7/8 has to be replaced by Pv($ = <^o) + Pv( < ^ ) = Vis))- The only 
thing that matters is the support of $ defined by 

Supp($) = {j : P A r($ = <pj) > 0} = {j : P*($ = ipj) > 0}. 

In this section we will generalize the concept of a wall. An easy adaptation 
of its proof will permit to draw the conclusions of Proposition [2] for a large 
majority of the 2 16 different supports. 

Consider an inhomogeneous BCA. We call adjacent cells {i + 1, . . . , i + p} 
an impermeable block of size p, if there exist b = (pi, . . . , b p ) G {0, 1} P , and 
(j i, . . . , j p ) G {0, . . . , 15} p such that if 

(x m (0), . . . ,x i+p (0)) = (6i, ...,b p ) and (cp i+1 , . . .,(p i+p ) = (tp Jl} . ..<pj r ), 
then for all t 

(Xi+x(t), . . . , Xi+p(t)), 

is unchanged, whatever values are chosen for Xi(s) and Xi +p+ i(s), for s = 
0, 1 . . . , t — 1. Here we assume that N is larger than p + 2. 
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As an example take p = 2, b = (0, 0) and (ji, j 2 ) = (2, 4). 
If a particular $ has (p^s in its support that yield an impermeable block, 
then any <3> having this set of ifj 's in its support will also have an imperme- 
able block. We therefore look for minimal sets of j's such that the tp^'s give 
rise to impermeable blocks. Here is a listing of these sets: 

g = {{0}, {1}, {7}, {8}, {14}, {15}, {2, 4}, {2, 5}, {2, 9}, {2, 12}, {2, 13}, 
{3, 4}, {3, 5}, {3, 10}, {3, 13}, {4, 9}, {4, 10}, {4, 11}, {5, 11}, 
{5, 12}, {6, 11}, {6, 13}, {10, 12}, {10, 13}, {11, 12}, {11, 13}, 

{2, 6, 10}, {4, 6, 12}, {9, 10, 11}, {9, 12, 13}}. 

Because of the symmetries of the collection of sets of CA rules (cf. Section[9]), 
of these 30 subsets there is a subcollection Q of only 12 where the dynamics 
of the associated is essentially different. Table [2] gives these supports 
with their ^-blocks ((p^ t , ■ ■ ■ y?j p ), and their 6-blocks (b\, . . . ,b p ). 

Table 2. Impermeable blocks. 



Support 


<^-block 


6-block 


Support 


<^-block 


6-block 


{0} 


(0) 


(0) 


{2,12} 


(2,12) 


(0,0) 


{1} 


(1,1,1) 


(0,1,0) 


{2,13} 


(13, 2) 


(1,0) 


{8} 


(8,8) 


(0,0) 


{3,5} 


(5,3) 


(0,0) 


{2,4} 


(2,4) 


(0,0) 


{3,10} 


(10,3) 


(0,0) 


{2,5} 


(5,2) 


(1,0) 


{10,12} 


(10,12) 


(0,0) 


{2,9} 


(2,9,9,2) 


(0,0,1,0) 


{2,6,10} 


(10,6,2) 


(1,1,0) 



Often there is more than one 6-block that will work, in the table we have 
chosen 6-blocks with minimal length. It can be quickly verified that the en- 
tries in Table[2]yield impermeable blocks. See Appendix [10] for two examples 
of the type of verification one has to perform. 

We can describe the supports which do not give rise to impermeable blocks 
by the list 

B = {{2, 3, 6}, {2, 3, 11}, {2, 10, 11}, {3, 9, 11}, {4, 5, 6}, {4, 5, 13}, 

{4, 12, 13}, {5, 9, 13}, {3, 6, 9, 12}, {5, 6, 9, 10}}. 

Any support which does not admit an impermeable block is a subset of a set 
in B. For an example of a proof that with these supports no impermeable 
blocks occur, see Appendix [TTJ 
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One can check that any subset of {0, 1, . . . , 15} is either a subset of a set 
from B, or contains a set from Q, and so we have completely described the 
supports with regard to the existence of impermeable blocks. 

Example. In [2] the two-connected RBN with solely rule ifr is studied, 
and it is claimed that "every node oscillates with period 2". Since if 7 is 
mirror equivalent to </?i, we see that actually the RBCA with support only 
consisting of <p 7 , has absorbing block (1,0,1). Let i be a node in the yj 7 -RBN 
with parents j and k. Then, if i happens to be one of the parents of both 
j and k (figure 8 configuration), and the state of i is 0, and that of j and 
k is 1, then these three nodes keep that state forever. This gives that the 
expected number of stable nodes in the RBN is at least 

N ■ F N (figure 8 configuration) • F N (initial values 1, 0, 1) = N 

5. Classifying RBCA's 

In the previous section we have established that all the RBCA's whose 
support contains a set in Q are well behaved in the sense that they are sta- 
bility continuous [o^ — > cr*), and do have a positive probability to stabilize 
(a* > 0). As there are only 53 supports described by B this means that this 
regular behaviour manifests itself for at least 99,9% of the RBCA's. The 
following result gives a complete classification. 

Theorem 1. An infinite one- dimensional random Boolean cellular automa- 
ton has no stable cells, i.e., cr* = 0, if and only if the support Supp($) of the 
rule distribution $ is a subset of{tp 2 , <Pto}, {fio, <fu}, {y? 4 , ¥12}, {<fi2, ¥13}, {^3, fe, ¥9, ^12} 
or {(f 5 ,ip 6 ,if 9 ,!f 10 }. 

Proof: From the results of the previous section it is clear that any CA with 
a* = must have its support in (the subsets of ) collection B. Because of 
symmetries it is enough to consider the 18 supports which occur as subsets 
of the sets in the collection B given by 

B = {{2, 3, 6}, {2, 3, 11}, {2, 10, 11}, {3, 6, 9, 12}}. 

In Section [7| we will prove that subsets of rules in {3,6,9, 12} all give sup- 
ports that belong to CA with a* = 0. 

It occurs that the remaining supports can be split into two parts. The first 
part consists of the supports which occur as subsets of B given by 

S = {{2, 3}, {2, 6}, {2, 11}, {2, 3, 6}, {2, 3, 11}, {2, 10, 11}}. 

The~ here indicates that we removed those supports that are symmetry 
equivalent to supports that are already listed, as, for example, {3, 11}. The 
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supports in S will be shown in Section [6] to be of the same regular type as 
found in Section HI The second part consists of the two supports 

{2} and {2,10}. 

Here the automata behave very similarly. First note that (p 10 is simply the 
left shift. The crucial fact is that 

(P2(x,y) = ip 10 (x,y) if (x, y) ^ (1, 1). 

First we consider the case Supp($) = {^p 2 }- Then ip 2 acts as the left shift 
for all times t > 1, since after one application of (p 2 there will be no blocks 
1*1 (i.e., 101 or 111) left anymore, since in the initial sequence 0*1 cannot 
overlap with itself with lag 2. The map (p 2 acts then as the shift on a subshift 
in which blocks of more than two ones can not occur. It follows that for each 
i the process (^Q(£)t>i) is the same subshift, which clearly can not stabilize. 
(Actually it can be shown that this process is a hidden Markov chain on 8 
equiprobable states.) 

In case Supp($) = {f 2 ,fio}, the CA runs as a left shift with occasional 
transformations of l's into 0's. This is a form of self-organization. In fact 
any block of l's that occurs at time to the right of a cell i with $j = ip 2 
will be cut into blocks of 1 and 11 separated by 0's. Since any cell i has 
a cell %' to the right of it with <Ev = f 2 with probability 1, there exists a 
(random) to, such that the process (Xi(t))t>t will only contain 1 and 11, 
and cannot stabilize. □ 




FIGURE 2. Evolution of the RBCA with support {^2} and {^2,^10} • 
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6. Stability by means of absorbing blocks 

Here we show that we may define a weaker notion than that of an imperme- 
able block which still permits us to conclude to continuity of stabilization 
and a positive probability to stabilize. 

Consider an inhomogeneous CA. We call adjacent cells {i + 1, . . . , % + p} 
an absorbing block of size p, if there exist b in {0, 1} P , and . . . , j p ) in 
{0, . . . , 15} p such that if 

(x m (0), . . .,x i+p (0)) = (6i, . . . ,b p ) and (<p i+1 , . . . , (p i+p ) = (<p 3l , . ..<p Jr ), 

then for some c with 1 < c < p the value Xi +C (t) is unchanged for all t, 
whatever values are chosen for Xi(s) and x i+p+ i(s), s — 0, 1 . . . ,t — 1. (Here 
we assume that N is larger than p + 2). The idea is that the evolution of 
one or both adjacent cells % and i +p + 1 may 'penetrate' the block, but will 
never influence the 'central' cell % + c. 

We give an example for an inhomogeneous automaton with rules y?2 and 
</?u. Let p = 4, b = (0,1,1,0), and (j 1; . . . , j 4 ) = (2,2,11,2). Then this 
block is absorbing with c = 2: 



cell % 


cell % + 1 


cell % + 2 


cell i + 3 


cell % + 4 


cell % + 5 




y?2 




<Pn 






X 





1 


1 





y 


x' 


1 — X 


1 








y' 


x" 


1-x' 








y' 


y" 


x'" 








1 


y" 


y'" 







1 


1 








Note that at time t — 4 the block 6 = (0, 1, 1, 0) reappears, so cell i + 2 has a 
periodic evolution with period 4, independently of the evolution of cell % and 
cell i + 5. This establishes that (2, 2, 11, 2) is an absorbing block. Similarly 
it can be verified that (2, 2, 3) is an absorbing block for the automata with 
rules y? 2 and <p 3 . 

For the support {2,6}, we present an absorbing block family with stable 
cells: take y?-block 

(2,2,2,6,6,6,2,2,2,6,6,6,2,2,2,2) 

and 6-block family 

(0,0,l,l,0,l,0,l,0,l,l,0,x,y,z, w ), 

where x,y,z and w can have any value as long as both xz and yw are not 
11. As a consequence of the </?2 rule, neither xz nor yw can be 11 at any 
later time. One can then check, in a case by case analysis for all remaining 
possibilities for xz and yw (three possibilities each), that the 6-block family 
keeps reappearing, and that the fourth and seventh cell are stable. We find 
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that if yw = 01 or yw = 10, it takes 8 time steps for another member of 
the 6-block family to reappear, for all other possibilities, this takes 4 time 
steps. 

It follows from these observations that the rules with supports in the col- 
lection S considered in the proof of Theorem [T] do not have stable cells. 
We can not yet conclude that a* > for supports {ip 2 , ^3} and {(p 2 , </?n} , 
since in both cases the central cell has a period 4 dynamics. However, we can 
find larger absorbing blocks with a stable central cell c: for {2, 3} take p = 9, 
c = 6, <^-block: (2, 2, 3, 2, 3, 2, 2, 2, 3) and 6-block: (0, 0, 1, 1, 0, 0, 0, 0, 1). For 
{2, 11} take p = 10, c = 6, <^-block: (2, 2, 11, 2, 11, 2, 2, 2, 11, 2) and 6-block: 
(0, 1, 1, 0, 0, 0, 0, 1, 1, 0). The existence of such absorbing blocks with a stable 
cell implies that a* > 0. 

7. Affine chaos 

Here we will analyze the RBCA's whose support is a subset of {3, 6, 9, 12}. 




FIGURE 3. Evolution of an RBCA with support {3, 6, 9, 12}. 

The crucial observation is that not only these rules are affine, but they can 
all be written as 

(p 3 (x,y) = x + f(y,(pj), 
where f(y, (p$) — 1, y, y +1, for j = 3, 6, 9, 12 respectively. It follows from 
this that we can write for some function f t with 4t — 2 arguments 

(3) Xi{t) = Xi_ t (0) + f t (xi- t+2 {0), • • • , x i+t (0), (fi-t+i, • • • , (pt+t-i)- 
Indeed, for t — 1 we have 

Xi(l) = Xj_i(0) + /i(x m (0), ifi), 
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where fi(y,<Pi) = f(y,fi) as above. For t = 2 we obtain, applying this 
equation or its shift three times 

Xi{2) = + 

= Xi_ 2 (0) + fi(xi(0), <£>i_i) + /i(xj(0) + /i(x t+2 (0), </?i+i), 

Defining 

/ 2 (^i(0), . . . , ar i+2 (0), v?i-i, • • • , yi+i) 

= /i(^i(0), + /i(a?i(0) + /i(x i+2 (0), v^i+i), 

we obtain Equation ([3]) for t = 2. Continuing in this fashion we obtain 
Equation (j3J) for all £. 

The following proposition will immediately imply that cr* = for the RBCA's 
from this section. 

Proposition 3. Consider an RBCA with Supp(<l>) C {y?3, (p^, <pg, ^12}, and 
let i be a fixed cell. Then the random variables pTj(t))^Q° are independent 
Ber(^) random variables under P*. 

Proof: We prove this by induction w.r.t. the length t of the cylinders. We 
will show that for each cell i and for all t — 1, 2, . . . 

(4) P*pQ(0) =e Q ,X i (l) =£!,..., Xi(t-1) =e t _!) = 2-* 

for all (e , £1, • • • , £t-i) from{0, 1}*. 

1 



z — t i — t + 1 z 2 + £ 

111 1 

AV t (0) £o 

\XXXXX)(XXXX/ 

xxxxx.xxxx/ 
xxxxxxx/ 

\ X X X X X 

\ / 

Xi{t) 



FIGURE 4. Additive propagation. 
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Although in general the (Xi(t)) will not be stationary (any cell which sta- 
bilizes but not at t — has a transient evolution!), a simple computation 
shows that Equation implies that shifted cylinders of length t also have 
probability 2~* to occur. 
We will use the following abbreviations: 

Y itt = f t {X^ t+2 (0), X i+t (0), . . . , $ i+t _i) 

E iM = {X i (p) = e ,X i (l) = e 1 ,...,X i (t-l) = e^ 1 }. 

First note that Equation PJ is true for t — 1 by definition of the RBCA. 
Figure H] is useful in the remainder of the proof. 

Suppose that Equation (J4]), i.e., F^E^g) = 2~* has been proved for cylinders 
of length t. Then for all Eq, . . . , s t -\ from{0, 1}* 

F.(E ittfe ,Xi(t) = 1) = 

P*(£ iit)£) JQ_ t (0) = 0,Y itt = 1) +P*(£ Mi£) X i _ t (0) = l,Y iit = 0) = 

p(x^(o) = o) p*(£ i)t , £ , y M = i) + p*(Xi_ t (o) = i) p*(£? M)6 , r i)t = o) = 
!(p*(£ iltl£ , = i) + p»(^ iAe , r <it = o) ) = 
|p,( J E; M , e ) = 2-(* +1 ), 

where we used that Xj_ t (0) is independent of £?i,t, e an d of 1^^. This equality 
proves Equation (TJJ for length t + 1 cylinders. □ 

What remains for the automata in this section is the question whether 
<tn —> For some automata this is easy to answer. We have already 
seen that for support {6} there is no stability continuity (see Proposition 
[1]). By symmetry, the same is true for support {9} and with a little more 
work this can also be shown for support {6,9}. On the other hand it is 
easy to show that &2N+i = and o~2N = 2 _iV for Supp($) = {3} and that 
<tat = 2~ N+1 for Supp($) = {12}, hence we do have stability continuity in 
these cases. Also for the combined support {3, 12} it is easily shown that 
°~n < 2~ N+1 , hence here too there is stability continuity. However, for the 
automata with supports {3,6}, {6,12}, {3,6,9}, {3,6,12}, {6,9,12}, and 
{3, 6, 9, 12}, the computation of ax becomes quite involved. We conjecture 
that there is stability continuity for these cases. 

8. Conclusion 

We have studied the behavior of one- dimensional random Boolean cellular 
automata with two inputs. Although this behavior can be quite diverse, 
we have shown that it does not depend on the probabilities with which 
the Boolean rules ip , . . . , (p i5 are attached to the integers, but only on the 
support of this random variable i.e. on the probabilities P($ = <^j) being 
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or positive. This contrasts with the behavior of random Boolean networks 
as determined by Lynch ([HI E]), see also [2J. His result is that there is 
ordered behavior as long as 

P($ = ip ) + P($ = <pui) > P($ = <Pe) + P($ = y 9 ) , 

and chaotic behavior when the opposite inequality holds. Nevertheless, there 
is an interesting parallel: for the one-dimensional random cellular automata 
P($ = (f ) + P($ = ipis) > implies regular behavior in the sense that cr* > 
and exponential convergence of cr/v to a*, while P($ = (/? 6 )+P($ = ^9) = 1 
implies that a* = 0, and no convergence of o~at to a*. 

9. Appendix: CA symmetries 

There are two symmetry operations on the collection of sets of CA rules. 
The first one is an extension of the mirror map Ai defined by 

M(0) = 1, M(1)=0. 

It is given by 

(M<p)(x,y)=M(<p(M(x),M(y))). 
The second one is space reversal TZ, defined by 

(n<p)(x,y) = <p(y,x). 

Both operations are involutions. The effect of Ai, written as a permutation 
is: 

(0 15)(1 7) (2 11)(4 13)(6 9)(8 14). 

The effect of 1Z is: 

(2 4) (3 5)(10 12)(11 13). 

10. Appendix: impermeable blocks 

Here we give two examples of the proofs that the blocks obtained from 
Table [2] yield impermeable blocks. 

To prove that the pair (61,..., 64) = (0, 0, 1, 0), (ji, . . . , j 4 ) = (2,9,9,2) 
yields an impermeable block, we have to consider a configuration as 





y?2 


<P9 


f9 


<P2 




X 








1 





y 



Filling in <f2{x, 0) = 0, <^g(0, 1) = 0, <^g(0, 0) = 1 and ^2(1, y) = 0, which is 
true for arbitrary x and y, we obtain the next line (time t — 1): 





¥2 


y?9 


y?9 


y?2 




X 








1 





y 


x' 








1 





y' 
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Since this is the same as the state at time 0, it follows that we do indeed 
have an impermeable block, consisting of four cells that always stabilize at 
time 0. 

There is only one exception of an impermeable block which does not consist 
of stable cells: 







y?3 




X 








y 


x' 


1 


1 


y' 


x" 








y" 


x'" 


1 


1 


y'" 



This time we obtain a period 2 impermeable block, but actually a period 1 
impermeable block is also possible by taking (bi, b 2 ) = (0, 1). 

11. Appendix: no impermeable blocks 

We will show that the RBCA's with their support in {{p 2 , f3, <Pn} do not ad- 
mit any impermeable blocks. We will do this by showing that diagrams as in 
the previous appendix can not exist. In the following we will use frequently 
that a block is impermeable if and only if its mirror image, respectively 
space reversal is impermeable. 

Consider the first element ip^ of the </?-block. This can not be ip 3 , since 
ip 3 (x,y) = x, which is not compatible with impermeability from the left. 
So it has to be tp 2 or tp n . Since -M(y?2) = ty?n, and M.((p 3 ) = ip 3 , we can 
assume without loss of generality that it is ip n . Since (pn(x,y) = if and 
only if (x, y) = (1, 0) the first element of the iterates of the 6-block must be 
equal to 1: 





<Pn 






V- 




X 

x' 
x" 


1 
1 










We next consider all three possibilities for ip i2 , filling in, if possible 
values of the second cell given by the rules. For j 2 = 2 we obtain: 




<Pu 










X 

x' 
x" 


1 
1 












However this gives a contradiction at t — 3, since <pn(x", 0) = 1—x" depends 
on x" . Exactly the same contradiction occurs when j 2 = 3. Consequently 
we must have j 2 = 11. The whole second column of the 6-block must consist 
of l's, again because of y?n(x, 0) = 1 — x. But then, since y?ii(l, 0) = also 
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the 3 column must be filled with l's starting from t — 1. It is then quickly 
checked in the following diagram that j 3 = 2 and j 3 = 3 are impossible: 





<pu 


¥>11 








X 




1 








x' 


1 


1 


1 






x" 


1 


1 


1 







Conclusion: also J3 = 11. Continuing in this fashion we find that the (string 
of indices of the) </>block has the form (ji, ■ ■ ■ , j p ) = (11,. ..,11). But then 
we have a problem in the p th column since (pn(l,y) = y depends on y. 
Hence for any p this possibility is ruled out. 
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